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IONOSPHERIC NON-LINEARITY 



SUMMARY 

The mechanism which causes the ionosphere to behave as a non-linear 
transmission medium when illuminated by a strong radio wave is discussed. 
It is shown that the non-linear effect is accentuated near the gyromagnetic 
frequency if the extraordinary wave is strongly-excited by the transmitter. 
Calculation of the magnitude of the non-linear effect by means of aray-'tracing 
technique is described and expressions are derived which enable the attenua- 
tion, distortion and cross-^modulation caused by non-linearity to be calculated. 



LIST OF PRINCIPAL SYMBOLS 



A 
Ao 



Go 



Eu 

e 
G 

M 
M' 
m 
N 
V 
Vo 



Total ionospheric attenuation of wave (dB) 
Total ionospheric attenuation of unmodula- 
ted wave (dB) 

Total ionospheric attenuation of wave in 
absence of disturbing wave (dB) 
Attenuation per unit path length of wave 
(dB/km) 

Attenuation per unit path length of wave in 
absence of disturbing wave 
Peak field strength of wave (mV/m) 
Plasma field strength, defined in Equation 
(5) (mV/m) 

Peak field strength of unmodulated wave 
(mV/m) 

Charge of electron 
A constant defined in Equation (2) 
Modulation index of wanted wave 
Modulation index of unwanted wave 
Mass of electron 

Number of electrons per unit volume 
Carrier amplitude of received signal 
Carrier amplitude of received signal in ab° 
sence of disturbing wave 



V 
Oo 

z 
dA{t) 



8A 



5A'u 

I 

F 



M' 



-A 

ay' 
m 



Average of random electron velocities 
Average of random electron velocities In ab- 
sence of disturbing wave 
Distance measured along path of wave (km) 
Instantaneous increase in attenuation due to 
ionospheric non-linearity 
Increase in attenuation of unmodulated wan- 
ted wave due to ionospheric non-linearity 
Increase in attenuation of unmodulated un- 
wanted wave due to ionospheric non-linearity 
Reduction of the fundamental component of 
the wanted wave modulation (dB) 
Reduction of the fundamental component of 
the unwanted wave modulation (dB) 
Electron collision frequency 
Electron collision frequency in the absence 
of disturbing wave 
tan-^UoyJGuo) 
Xan^(2GyJGuo) 

Angle between direction of wave and Earth's 
magnetic field 
Angular frequency of wave 
Angular gyromagnetic frequency 
Angular modulation frequency of wanted wave 
Angular modulation frequency of unwanted 
wave 



1. INTRODUCTION 

Effects caused by non-linearity in the ionos- 
phere are becoming increasingly important with the 
trend towards the use of higher powers for night- 
time sky-wave broadcasting on low and medium 
frequencies. These effects may be classified as 
follows: 



(11) the ionospheric absorption of a strong wanted 
wave may be increased and harmonic distor- 
tion of the audio-frequency modulation may 
occur. 

(iii) when two strong waves are radiated simultan- 
eously, intermodulation may take place in 
addition to cross-modulation. 



I) the modulation of a powerful disturbing wave 
may be impressed on another wave traversing 
the same region of the ionosphere, so caus- 
ing interference. This is the so-called Luxem- 
bourg effect, first observed in 1933. 



There is considerable evidence, both theoretical 
and experimental, which suggests that these effects 
may be enhanced if the frequency of the disturbing 
wave is close to the gyromagnetic frequency. 



The general theory of ionospheric non-linearity 
is surveyed in this report and extended to the 
situation which arises in tropical latitudes, where 
the extraordinary wave is excited strongly. Use is 
made of a simple ray-tracing technique to assess 
the effect of non-linearity on an incident wave 
which excites the ordinary and extraordinary waves 
in roughly equal proportions. The dependence of 
distortion and cross-modulation on the modulation 
frequency is discussed and the theory extended to 
include effects occurring when two strong waves 
illuminate the same region of the ionosphere. 



to each electron. Now the constituents of the ionos- 
phere are in continuous random motion and collisions 
between electrons and heavier particles occur fre- 
quently. When such a collision occurs the ordered 
energy supplied ty the wave is converted into 
random energy and more energy has to be supplied 
by the wave to set the electron in motion again; 
the wave therefore suffers loss and is attenuated. 
The more frequently the electrons collide the greater 
is the attenuation, unless collisions occur so often 
that the electrons have insufficient time to acquire 
the full amount of ordered energy which the wave 
can supply before the next collision takes place. 



2. THE MECHANISM OF IONOSPHERIC NON- 
LINEARITY 

The ionosphere consists of a mixture of neutral 
gas molecules, positive ions and free electrons. 
The electrons, being relatively light, are easily 
deflected by electric and magnetic fields. When a 
wave propagates through the ionosphere its electric 
field imparts an oscillatory motion to the electrons, 
which re-radiate like miniature aerials. The re- 
radiation modifies the. velocity of propagation of the 
wave and refraction occurs when the electron den- 
sity varies, so that if the frequency is not too high 
waves are refracted back towards the Earth. 

If the Earth's magnetic field were absent the 
oscillatory motion of the electrons would be parallel 
to the direction of the electric field of the incident 
wave and the re-radiation would have the same 
polarization as the incident wave; the wave would 
therefore traverse the ionosphere without change 
of polarization. The Earth's magnetic field modi- 
fies the oscillatory motion of the electrons, causing 
them to move in complicated orbits. Their re- 
radiation is not of the same polarization as the 
incident wave and consequently the polarization 
of the resultant wave changes continually as it 
traverses the ionosphere. 

There are, however, two incident-wave polar- 
izations (In general elliptical) which excite the 
electrons in such a way that the waves re-radiated 
by the electron have the same polarization; these 
waves therefore traverse the ionosphere without 
change of polarization. These two special polar- 
izations are termed the 'ordinary' and 'extraordinary' 
waves, the parameters describing their polarizations 
depending on the angle between the direction of 
propagation and the direction of the Earth's mag- 
netic field. Since these two waves traverse the 
ionosphere with different velocities and attenua- 
tions it is usual to resolve the incident wave into 
the two components and then consider them separ- 
ately. 

Both the ordinary and extraordinary waves 
impart a minute but significant quantity of energy 



Up to this point we have been considering 
waves which are so weak that the oscillatory motion 
which they impart to the electrons is small com- 
pared with the random motion. If the wave is suf- 
ficiently strong, however, the random motion of the 
electron may be increased significantly by transfer 
of energy from the oscillatory motion at each col- 
lision. The chances of a collision occurring are 
then greater, so that the collision frequency is 
increased and this may result in an increase of 
absorption, as discussed in the previous para- 
graph. This increased loss would apply not only to 
the strong wave but to all other waves traversing 
the same region of the ionosphere, however weak 
they might be. If the strong wave (which wilt be 
referred to as the disturbing wave) is amplitude- 
modulated at a sufficiently low audio frequency the 
attenuation will vary at audio frequency and the 
modulation of the disturbing wave will be super- 
imposed on the weaker wave; this is the basis of 
the Luxembourg effect. If the disturbing wave itself 
is being received, a non-linear relationship will 
exist between the received and transmitted powers 
and its waveform will be distorted as a result. 

Special consideration must be^ given to the 
situation which arises when the disturbing wave 
excites the extraordinary wave. Near the gyro- 
magnetic frequency a relatively weak extraordinary 
wave can impart a large ordered motion to the elec- 
trons and cause a significant increase in collision 
frequency. This gyro-resonance effect was proposed 
theoretically by Bailey^'^ in 1937 and has since 
been demonstrated experimentally-^''^. Here we 
consider the situation which arises if the disturb- 
ing wave is the extraordinary component of a wave 
whose frequency is near the. gyromagnetic frequency, 
the wanted wave being the ordinary component. 
Ordinary and extraordinary waves are excited in 
roughly equal proportions on East-West paths at 
magnetic dip latitudes of about 45° and on long 
North-South paths near the magnetic equator. 

To determine the influence of a disturbing 
wave on a wanted wave, the increase in collision 
frequency due to the disturbing wave is first calcu- 
lated. The change in the absorption coefficient of 



the wanted wave is then found; this yields a factor 
by which the attenuation of the wanted wave must 
must be multiplied. Provided the attenuation of the 
wanted wave in the absence of the disturbing wave 
is known, the change in attenuation can be found. 
It is then possible to calculate the effect which 
the attenuation change has on the modulation of 
the wanted wave; the cross-modulation impressed 
on other waves may also be determined. 

The theory of ionospheric* non-linearity is 
discussed in more detail in the sections which 
follow. 



When a collision occurs this energy is con- 
verted into random energy and is shared between 
the electron and the heavy particle with which it 
collides. The mean random velocity of the elec- 
trons therefore increases from vq In the absence of 
the wave to u in the presence of the wave. When a 
steady state has been reached the average energy 
of an electron is increased by an amount given by 



V2 m(v'' - vg)G 



w 



2) 



where G is a constant which has been determined 
experimentally** and is approximately equal to 
1.7 X 10"^. 



3. INCREASE IN COLLISION FREQUENCY DUE 
TO A DISTURBING WAVE 

3.1. Increase in Collision Frequency when the 
Earth's Magnetic Field is Neglected 

Since the constituents of the ionosphere are 
in continuous random motion, collisions frequently 
occur. If we follow the movement of a single elec- 
tron we find that, on average, v collisions occur 
every second even though the intervals between 
successive collisions are irregular.* In discuss- 
ing the behaviour of the ionosphere the effective 
collision frequency is assumed to be slightly greater 
than the actual collision frequency, for reasons 
which are beyond the scope of this report. In all 
that follows, 1/ will be assumed to denote the effec- 
tive collision frequency. 

If the Earth's magnetic field is neglected it 
may be shown that the work done between collisions 
by the wave in setting an electron in motion is, on 
average, given by 



1 e' 

10= - — 



E' 



(1) 



per collision, where 



e and m are the charge and mass of the elec- 
tron, respectively 

A 

E is the peak value of the incident electric 
field 

u is the collision frequency 

oj is the angular wave frequency 



* In the D- and E-regions, where most absorption takes 
place, collisions with neutral molecules occur much 
more frequently than collisions with positive ions. If on 
electron collides with a positive ion its free existence 
may be terminated and it then plays no further part in the 
propagation of the wave. The concept of on average 
collision frequency is, however, still valid for the free 
electrons which remain. 



Now the collision frequency u is directly 
proportional to v because the mean free path travers- 
ed by the electrons between collisions is independ- 
ent of their velocity. From Equations (1) and (2) 
we may therefore obtain the following result 



W Vv^^'W'^o^ "'" 



(3) 



to 



where uo is the collision frequency in the absence 
of the incident wave. Equation (3) shows that u 
increases in the presence of the wave. 



Equation (3) may be written in the form 



where 



Ep = vo 



^ G(vo +c^ )\ 



(4) 



(5) 



The quantity Ep has the dimensions of field strength 
and is known as the "plasma field strength*; in the 
D- and E-layers it is given approximately by 



E. 



150 




V2 



mV/m 



(6) 



where 1^0/277 and /"are in MHz. At heights greater 
than 80 km, tL>^»7y^ throughout the m.f. band; the 
expression for E t'^en simplifies to 



EpCil50/ mV/m 

In the m.f. band, E/Ep is unlikely to exceed 
unity even with the highest powers contemplated. 



**There is some uncertainty about the correct value for 
and its value cannot be derived from simple theoretical 
considerations. For a discussion of this problem see 
Appendix II of Reference 4. 



Consequently the increase in v is small compared 
with CO and Equation (4) may be simplified to 



4 



1 + (E/Ep) 



with little error; this approximation was made by 
Huxley and Ratcliffe.^ In Section 3.2 we shall 
show that a modified form of Equation (4) applies 
when the Earth's magnetic field is taken into con- 
sideration and that the quantity equivalent to E/Ep 
may be appreciably greater than 1. The approximate 
solution of Huxley and Ratcliffe is no longer valid 
and we must use the exact solution of Equation (4), 
which is a quadratic in i^^. For convenience the 
exact solution will be stated here; it is 




2/2 2\ 

CO + [V^ + CO ) 






(7) 



This result has been given in slightly different form 
by Ginzburg and Gurevich.*^ Curves showing the 
relationship between v/i^o and E/Ep derived from 
Equation (7) are shown in Fig. 1. 



wave, M the modulation index and co^ the angular 
modulation frequency. When Huxley and Ratcliffe's 
approximation is valid the ionosphere may be repre- 
sented by the electrical analogue model illustrated 
in Fig. 2. This shows that v^ rises and falls with 
time constant 1/Gi^o, in the E-layer approximately 
equal to 1'25 milliseconds. At very low modulation 
frequencies, v varies over the modulation Cfc\e in 
accordance with the *steady-state' solution. At 
high modulation frequencies, however, only the 
constant term in^he expression for E^ is important; 
this is equal to Eu(l + M^/2), which is proportional 
to the total power contained in the carrier and side- 
bands. Thus at high audio frequencies v remains 
constant over the modulation cycle, its value depend- 
ing only on the power density of the incident wave. 
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Fig, 1 " Increase of collision frequency as a func- 
tion of strength of disturbing wave 

Earth's magnetic field neglected 

It should be noted that the values of i^/i^o given 
by Equations (3) and (7) are 'steady-state* values, 
applicable when E is constant, i.e. when the dis- 
turbing wave is ur^modulated. If the wave is ampli- 
tude modulated E varies with time and may be 
written 

A A 

E = Eu (1 + !Vf cosoj^O 
•where Ey is the peak value of the unmodulated 



Fig. 2 - Analogue model of the ionosphere 

Time constant CR = 1/Gl^o 



Input voltage^ £(/) 
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+ 2McosCO^t +— cos 2CiJ^t 



The effect which the variation of u at audio 
frequency has on the received signal is discussed 
fully In Section 6. 

3.2. Increase in Collision Frequency, taking 
into account the Effect of the Earth's Mag- 
netic Field 

When the effect of the Earth's magnetic field 
is taken into consideration the average energy 
acquired by the electrons between collisions is 
modified and Equation (1) applies only when the 
electric vector of the incident wave is parallel with 
the direction of the Earth's magnetic field. The 
theory contained in Section 3.1 must therefore be 
restated with a modified value of w. 

2 7 

It has been shown by Bailey that the average 
work done by a plane-polarized incident wave on an 
electron is 



w = 



4m 



I L^ + {co~coy^f 

^ 1 1 . . 

+ ™ — 2 Sin 



(8) 



-A 



2 cos lA 



I/^ +CO- 



where oj^ is the angular gyromagnetic frequency 
and ip is the angle between the directions of the 
electric field of the incident wave and the Earth's 
magnetic field. We may compare this equation with 
Equation (1); as expected, Equation (8) reduces to 
Equation (1) when i// = 0, since the Earth's mag- 
netic field then has no effect on the motion of the 
electrons. When i/^ f 0, the first term of Equation 
(8) is large at frequencies near the gyromagnetic 
frequency, because of the resonant motion of the 
electrons. 

We now proceed as in Section 3.1, replacing 
Equation (1) by Equation (8). The complete expres- 
sion for ly/i/Q leads to a quartic in i/ , Bailey has 
obtained a solution valid for relatively weak dis- 
turbing waves by replacing i/ in Equation (8) by vq] 
this is equivalent to the approximation of Huxley 
and Ratcliffe. There is, however, less justification 
for making the approximation in this case because 
large increases in v can occur near the gyromag- 
netic frequency. 

A new approach to the problem is adopted 
here by restricting the calculation of 1//1/0 to ^/^ = 
and ip = 90° only. 

When v^ = the result given in Section 3.1 
applies. When i/j = 90° the last term on the right- 
hand side of Equation (8) is zero and it can be 
shown that, in the m.f. band, the second term in the 
square bracket is always small compared with the 
first; the second term will therefore be neglected. 

Substitution of the resulting value of w in Equation 
(2) then leads to the following quadratic in 1/^, 
equivalent to Equation (3) 



still be used to estimate how u varies with modula- 
tion frequency provided E^ is replaced by Ep. 

In the D-layer there is little difference between 
the values of u/uo for i^ = and 90° and therefore 
no great uncertainty about the value of ly/i^o for 
other values of i//. In the E-layer it is convenient 
to calculate 1//1/0 for ordinary and extraordinary 
disturbing waves separately; for transverse propaga- 
tion the two waves correspond to if/ ^ and i// = 90° 
respectively. The value of 1//1/0 due to the com- 
bined effect of the two disturbing waves cannot 
easily be calculated; the larger of the values for 
the individual waves is used here as a basis for 
further calculation. 



4. CHANGE IN ABSORPTION OF WANTED WAVE 
DUE TO INCREASE OF COLLISION FRE- 
QUENCY 

The relationship between absorption and col- 
lision frequency is complicated and has been dis- 
cussed fully in Chapters 9 and 10 of Reference 8. 
At medium frequencies, however, certain simplify- 
ing assumptions may be made for paths involving 
oblique incidence, since most absorption occurs at 
levels where the refractive index does not differ 
greatly from unity. Here the attenuation of the 
ordinary wave is proportional to i//(iy^ + co'^) provid- 
ed the direction of propagation is perpendicular to 
the direction of the Earth's magnetic field. The 
rate of attenuation a of a wanted wave is then 
related to its rate of attenuation oq in the absence 
of the disturbing wave by the expression 



(t; ^''2y\,^G.% 



(A^) 



(9) 



a ly ly^ + CO iv 1 + (cl)/i^o) 



Oo Md 1/2 + &j2 i/Q [v/vo)'^ + (w/md)2 



(12) 



where Aoj - |6j-cl),_,|, the deviation from the gyro- 
magnetic frequency. 



We now define a 'modified plasma field-strength* 



Epas 



e'p = v,~^g[i . (^o.)^'' = E^ m^s^^ 



\y-2 



Equation (9) may then be written as follows 



2 2 



2 



v^ + (A&>)2 



(IJ 



(10) 



11) 



Comparison of Equations (4) and (11) shows that 
they are of the same form. Consequently vIvq may 
be determined from Fig. 1 if Ep is replaced by Ep 
and (o by Acu, The equivalent circuit of Fig. 2 may 



Thus if v/vq is known, aloo can be determined. 

Fig. 3 shows how aloo varies as vIvq increas- 
es. It will be seen that if vo is greater than oj the 
attenuation always decreases, for reasons which 
were discussed in Section 2. If vq is less than co 
the attenuation increases until i/ = co, after which it 
falls. 



5. THE NIGHT-TIME BEHAVIOUR OF THE ION- 
OSPHERE AT MEDIUM FREQUENCIES 

Although the theory discussed in the previous 
sections enables the fractional increase in attenua- 
tion caused by disturbing wave to be calculated, 
it is of no practical use unless the attenuation 
which occurs in the absence of the disturbing wave 
is known. To assess this, the properties of the 
ionosphere must be studied more closely. 
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Fig. 5 - Change of attenuation of wanted wave as 
a function of increase in collision frequency 

Transverse propagation of the ordinary wave 



Figs. 4(a) and (5) show typical night-time 
electron density and collision frequency profiles 
for the D- and E-layers, taken from Figs. 1 and 7, 
curve (5), of Reference 9 but with the maximum 
electron density increased to 2000 per cm^ as this 
value is believed to be more accurate. ^^ Fig, 4(c) 
shows the refractive index for transverse propaga- 
tion of the ordinary Wave at 1 MHz. It will be seen 
that the refractive index is very close to unity in 
the Olayer and does not fall below 0-915 in the 
E-layer. Thus 1 MHz waves incident at angles less 
than arc sin 0-915 (66^ will penetrate the E-layer 
at night; the skip distance is 440km. 
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(a) Electron density (b) Effective collision frequency 

(c) Refractive index for transverse propagation of the ordinary wave at 1 MHz 

(d) Vertical-incidence attenuation for transverse propagation of the ordinary wave at 1 MHz (one traversal of the layer) 



The rate of attenuation of a vertically-incident 
1 MHz ordinary wave is shown in Fig. 4(d). It is 
approximately proportional to Nv where N Is the 
electron density, and rises to a maximum, as N 
increases and v decreases, at about 74km. This 
occurs In the Olayer, which at night should be 
regarded more as a region of enhanced absorption 
than as a region of increased electron density. The 
sudden rise in electron density in the E-layer 
causes the rate of attenuation to rise to a second 
and greater maximum at about 94 km, despite the 
greatly reduced collision frequency at this height. 

In calculating the effect of a disturbing wave 
It is important to distinguish between D- and E- 
layer attenuation, because vo differs by a factor 
of more than 20 between the two layers. The calcu- 
lation of D-layerloss presents no difficulty because 
at m.f. the wave is not deviated significantly and 
the path may be assumed to be a straight line 
through the layer, even at oblique incidence. The 
E-layer loss is not so easy to determine because it 
depends on the path length and depth of penetration 
in the layer. A crude estimate may, however, be 
made from propagation curves based on measured 
field strengths and published by the EBU^^ by 
subtracting the calculated D-layer loss from the 
total ionospheric attenuation. For more accurate 
estimates use must be made of ray-tracing tech- 
niques, as described in Section 7. 

Fig. 5 shows D-layer losses calculated from 
the profiles of Fig. 4 and E-layer losses estimated 
from the EBU curves, and shows how these losses 
vary with frequency and distance. They represent 
median values for night time. Although the D-layer 
attenuation decreases as the frequency rises the 
E-layer losses show the opposite trend; despite the 
decreased rate of attenuation at the higher frequen- 
cies, waves penetrate more deeply into the E-layer 



and consequently suffer a greater total attenuation. 

5.1. Non-linearity in the D-Layer 

The absorption of the wanted wave is greatest 
at 74km, where the collision frequency vq is 3-7 x 
10^ per sec. Thus t/q = cd at about 0-6 MHz. Fig. 3 
shows that an increase in v will cause the absorp- 
tion to decrease at frequencies below 0-6 MHz and 
suggests that a/oo is unlikely to exceed 1*3 at any 
frequency in the m.f. band. The highest values of 
a/oo occur at the highest frequencies, but here the 
total D-layer loss is least, as Fig. 5 shows. We 
therefore conclude that non-linearity effects are 
unimportant in the D-layer; this conclusion has 
been confirmed by detai led calculation for an 800 km 
path. 

5.2. Non-linearity in the E-Layer 

On entering the D-layer the incident wave may 
be resolved into ordinary and extraordinary waves 
which suffer different attenuations. When the 
ordinary wave reaches the E-layer it may still be 
sufficiently strong to act as a disturbing wave. 
The extraordinary wave may also act as a disturb- 
ing wave even though it is more heavily attenuated 
in the D-layer because of gyro-frequency resonance. 
Here the effect of each wave is considered separate- 
ly because the theory for assessing their combined 
effects has not yet been developed. 

In the E-layer at m.f. oj»vo and consequently 
a/oo may be assumed to be equal to v/vo\ this 
result follows from Equation (12). 

5.3. Non-linearity on a Typical Path 

Although non-linearity is important only in the 
E-layer, allowance must be made for attenuation in 
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Fig. 5 - Variation of ionospheric attenuation with distance and frequency 
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the Olayer. Fig. 6 shows verticahincldence losses 
for ordinary and extraordinary waves, calcuTated by 
integrating curves similar to Fig. 4{d). If these 
josses are divided by the cosir>e of the angle of 
incidence, the strengths of both waves at the base 
of the E-layer may be calculated. Values of i//vo 
and hence o/oo, may then be calculated for the base 
of the E-layer, using the method described in Sec- 
tion 3. 



3 5 
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fraquency, MHz 

Fig, 6 "Vertical-'incidence attenuation in the Delayer 

Transverse propagation. Gyromagnetic frequency 920 kHz, 



Fig. 7 shows the values of a/oo which would 
occur at the base of the E-layer at distances of 
400 and 800km from transmitters having effective 
radiated powers (e.r.p.s) of TVa and 30 MW. The un- 
broken curves are the values of a/oo occurring when 
the transmitters excite only the transverse ordinary 
wave; these curves apply to propagation in ail 
directions near the magnetic poles and are indepen- 
dent of the gyromagnetic frequency. The unbroken 
curves are approximately correct for European 
latitudes since the excitation of the extraordinary 
wave by a vertical aerial is relatively weak. 

The broken curves show the values of a/oo 
occurring if the transmitters excite the transverse 
extraordinary wave strongly; this situation would 
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- Increase in attenuation of ordinary wave 
at base of E4ayer 

(a) Reflection point 400km from transmitter 

(b) Reflection point 800km from transmitter 



— __ Transverse ordinary di sturbing wav@ 
— — — — Transverse extraordinary disturbing wave; 
gyromagnetic frequency 920kHz 



arise when transmitting on an East-West path near 
the magnetic equator. An interesting feature of the 
broken curves is the double-hump shape centred on 
the gyromagnetic frequency, the central minimum 
being due to the high absorption loss suffered by 
the wave in passing through the D-layer. It is of 
interest to note that double-hump curves with 
roughly the same frequency separation between 
maxima have been observed in cross-modulation 
experiments by Cutolo in Europe-^ and Hibberd In 
Australia.'^ 

The large values of a/oo due to the extra- 
ordinary wave apply only at the point where the 
wanted wave enters and leaves the E-layer, because 
the extraordinary wave is rapidly attenuated above 
the base of the layer. Furthermore, for frequencies 
above the gyromagnetic frequency, the extraordinary 
wave is reflected at a lower level than the ordinary 
wave and never reaches the upper part of the ordin- 
ary-wave path, where most of the ordinary-wave 
attenuation occurs; in this region, therefore, any 
non-linear effects are caused by the ordinary wave 
alone. 



The overall increase in the attenuation of the 
ordinary wave radiated from a high-power trans- 
mitter must therefore be calculated by tracing the 
path of the wave through the layer and finding the 
attenuation increase at each point. This may be 
done by regarding the E-layer as a set of uniform 
horizontal slabs each 1 km thick and using SnelTs 
Law to trace the ray paths. As an example, Fig. 8 
shows paths calculated in this way, drawn as 
smooth curves with an exaggerated vertical scale. 
The Appleton-Hartree formula^ was used to calcu- 
late the refractive index and absorption index in 
each layer, transverse propagation being assumed 
throughout. 

The ordinary and extraordinary waves whose 
paths are shown in Fig. 8 were assumed to be ex- 
cited in roughly equal proportions by the transmitter; 
this situation would arise, for example, near 45° 
dip latitude on East-West paths. The attenuation 
of the ordinary wave, when a low-power transmitter 
is used, was calculated by adding attenuations 
calculated from the Appleton-Hartree formula for 
individual layers. The additional attenuation 
resulting from the use of a high-power transmitter 
was then estimated as follows. The strengths of 
the ordinary and extraordinary interfering waves in 
each layer were calculated and the values of vIvo 
and hence a/ ao, which each would cause individual- 
ly were determined from Fig. 1; the larger of the 
two values was adopted where both waves are 
present because of the difficulty, mentioned earlier, 
of calculating the combined effect of the two waves. 
The weak-signal attenuations previously calculated 
for the ordinary wave for each layer were multiplied 
by the appropriate values of al Oo and again added, 
so giving the total increase in attenuation due to 
non-linearity. It is worth noting that the attenuation 
increase suffered by the wave on its downward path 
is similar to the increase on the upward path even 
though the wave has been attenuated in the E-layer, 
because other upgoing waves, shown in Fig. 8, dis- 
turb the downcoming wave. 

Calculations were made for an 800km path at 
frequencies of 1.1, 1-3 and 1-5 MHz. The calcula- 
ted weak-signal E-layer attenuations at these fre- 
quencies are 6-9, 7-7 and 13-1 dB respectively; 
these losses are consistent with those estimated 
from the EBU curves, shown in Fig. 5(b), the high 
value at 1-5 MHz arising because the ordinary wave 
almost penetrates the layer at this frequency. The 
attenuation increases due to non-linearity were 
calculated for three values of e.r.p. and were found 
to be similar on all three frequencies. Fig. 9 shows 
average values for the three frequencies plotted as 
a function of e.r.p.; it will be seen that the rela- 
tionship between the e.r.p. and the attenuation 
increase is approximately linear for the range of 
powers considered. On theoretical grounds a linear 
relationship would not be expected unless Huxley 
and Ratcliffe's approximation were valid. The fact 
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Fig. 9 - Increase of attenuation due to non-'linearity 
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40° dip latitude; 800 km path; Frequency M to 1.5MHz; 
Gyromagnetic frequency 920 kHz 

that the calculation described in this section leads 
to a linear relationship shows that Huxley and 
Ratcliffe's approximation is justified for the path 
as a whole even though It may not be valid where 
the wave enters and leaves the E-layer. The cal- 
culations showed that the greater part of the ord- 
inary-wave attenuation occurs near the apex of the 
trajectory and here the approximation is valid for 
the powers considered. 



6. AUDIO-FREQUENCY DISTORTION 
BY IONOSPHERIC NON-LINEARITY 



CAUSED 



In this section we consider how the modulation 
of the wanted wave is affected by the disturbing 
wave. Huxley and Ratcliffe's approximation will 
be used, since it has been shown In Section 5 to 
be valid for a typical path taken as a whole. 

Suppose that z represents distance measured 
along the path of the wanted wave and that ao{z) is 
its rate of attenuation, in dB per unit path length, 
in the absence of the disturbing wave. The total 
attenuation is therefore Ao = fao{z)(^z. 

When the disturbing wave is present the rate 
of attenuation increases to a{z). We may therefore 
write 

^ -■ (13) 

provided E{z) varies very slowly with time. In 
Equation (13), E(z) Is the peak value of the disturb- 
ing wave and Ep{z) the plasma field strength, or the 
modified plasma field strength in the case of an 
extraordinary disturbing wave. Both quantities 
vary with height and are therefore functions of z. 



10 



If E(z) varies more rapidly with time, the in- 
stantaneous rate of attenuation depends also on the 
ionospheric time constant I/Guq. If 

E(z) =^u(^) ^^ +^ oosco^t] 

then the electrical analogue mode! of Fig. 2 shows 
that i^^ is given by 



2 2 



fr) 



at) 



where 



fit) =1 +^ +2MCOS01 cos(a>^^-0i) + 



+ ™ COS02COS(2c^rn^-02) (14) 

<i>i = ^ar^' Hoj^CR) ^ tan' \co J Gi^o) 

<h2 = tan" ^{2oj^CR) ^ tan" ^{2coJGuo) 

In the E-layer the attenuation increase is therefore 
given by 

a{z) A^) 



ooiz) ^(^(z) 




(15) 



Integration of Equation (15) gives the total 
instantaneous attenuation of the wanted wave, as 
follows: 



^(0 



-Aa(z)6z = (a,,(z)6z + Mt) ( oo (z) (f^j 



6z 



(16) 



The first term on the right hand side of Equation 
(16) is the attenuation which would occur in the 
absence of the disturbing wave and the second 
term is the increase due to non-linearity.* Since 
f(t) = 1 when the disturbing wave is unmodulated, 
the second term may be written in the form 

8A{t) = f(t)8A^ 

where 8A(t) isthe instantaneous increase in attenua- 

The linear relationship between the attenuation increase 
and the e.r.p., referred to in Section 5, moy be demon- 
strated by making the substitution E\j{z) = a(z)Ey(0) 
in Equation (16)^ where a{z) is a function which de- 
scribes how B\j{z) varies along the path of the wanted 
wave. The attenuation increase is then seen to be 
proportional to the square of £^(0), the field strength 
at. the point where the wave enters the ionosphere; 
[E(j(0)] is proportional to the e.r.p. 



tion due to non-linearity and My is the increase 
which occurs when the disturbing wave is unmodula- 
ted. The calculation of 8A^ was discussed in 
Section 5. 

6.1. Cross-modulation 

The wanted wave is any wave, however weak, 
traversing the same region of the ionosphere as 
the disturbing wave. The variation in the attenua- 
tion of the wanted wave, due to the disturbing wave, 
causes the modulation of the latter to be super- 
imposed on the wanted wave. The strength of the 
wanted wave is also reduced. 

Let Vq be the carrier ampi itude of the received 
signal in the absence of the disturbing wave. When 
the disturbing wave is present Vq is decreased by 
8A{t) decibels to a new amplitude V(t), We therefore 
have 

20 log 10 i~ = 8-68 log ~ ^ 8A{t) 
Provided 8A(t) is less than 2dB we may write 



V(t) ^Vo exp - 



Vo 




^Vo 



1- 



8A(t) 
8.68 



where 8A\j is the increase in the attenuation of the 
wanted wave caused by the disturbing wave when 



the latter is unmodulated. 



*« 



Substitution of an expression for f(t)^ similar 
to Equation (14) then leads to the following results: 

(i) The strength of the wanted wave is reduced by 

SaV . m'u (l . f^)clB (18) 

where M is the modulation index of the dis- 
turbing wave. 



(ii) The modulation of the disturbing wave is 
superimposed on the wanted transmission with 
a level given by 

L*i = 20 logio(M" coS(/)i' My) + 5A'p -12-8 dB 

(19) 

(iii) The second harmonic of the disturbing wave 
is also superimposed with level 



|(MVcos02<5A'uJ 



l's = 20 logio((MVcos02M| 



+ 8Af -24.8 dB 
(20) 



Certain quantities in this section are primed to make 
the notation consistent with that of Section 6.3. 
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In Equations (19) and (20) the reference level 
(OdB) corresponds to 100% cross-modulation. 

6.2. Self-distortion 

The wanted wave and the disturbing wave are 
one and the same. At the modulation peaks the 
increase in attenuation due to non-linearity is 
greater than it is when the wave is unmodulated; 
at the modulation troughs it is less. Non-linearity 
therefore reduces the peak-to-trough amplitude of 
the received signal (assumed to be derived from a 
detector in the absence of automatic gain control) 
and also generates harmonic distortion.* 

In the absence of non-linearity the received 
signal would be Vo(1 + Mcoscoj^t); because of non- 
linearity it is V(t). We therefore have 



20 logio 



Vo(1 + fWcosco^O 



V{t) 



Hence 



y{t)^Vo{^ +Mcosco^O 



8A(t) 



8.68 



(21) 



provided dA(t) is less than 2dB. Now dA(t) is 
equal to f(t) SAy, where f(t) is given by Equation 
(14). If this value of 8A(t) is substituted in Equa- 
tion (21) and the latter expanded, the following 
approximate results are obtained, provided 8A\j 
does not greatly exceed 0*5 dB. 

(i) The carrier amplitude of the received signal is 
reduced by 



8A 






M^cos^0i 



dB 



(ii) The fundamental component of the received 
signal is reduced by 



8Af =8A^ 



1 + 2cos'^<?!)i +7 (2 + cos ^0^2) 



dB 



(22) 



(iii) The second-harmonic level (relative to the 
fundamental) is 



(iv) The third-harmonic level is 

La =. 20 togio(fV[^cos02 SAy) + 8Ap -30.8 dB 

(24) 

At very low modulation frequencies the fundamental 
is reduced by 3(1 + fW^/4) SAy dB; this is more 
than three times the reduction of the unmodulated 
wave, 8A\j» At high modulation frequencies, how- 
ever, the reduction is only slightly greater than 
8A\j, The reduction lies about midway between 
the two limiting values when co^ = Gvo i.e. when 
the modulation frequency is about 130 Hz. The 
harmonic distortion is greatest at low modulation 
frequencies and decreases at 6dB per octave at 
high modulation frequencies. 

The most important effect is the reduction of 
the fundamental component of the received signal 
because this is what determines the reduction in 
receiver signal-to-noise ratio. It is mainly due to 
the increase in the attenuation of the signal as a 
whole rather than to a reduction of the modulation 
depth. At high modulation frequencies the carrier 
and sidebands are reduced by exactly the same 
amount and there is no change in the modulation 
depth. Only at low modulation frequencies is the 
modulation depth reduced, but this reduction is 
small compared with reduction of the carrier fre- 
quency. Recent experiments^^ with high power 
l.f. transmitters have confirmed that demodulation 
is unimportant, as no significant difference be- 
tween the percentage modulation of transmitted and 
received waves was observed. 



6.3. Interaction between Two Strong Waves 

We consider now the situation which arises 
when the wanted wave and the disturbing wave are 
of comparable strength; this could occur, for exam- 
ple, when a station transmits simultaneously on 
two frequencies. 

We make the initial assunr.ption that the work 
done on an electron by two waves of different fre- 
quency is equal to the sum of the work which each 



L2 = 2Ologio{M[cos^0i + cos0i cos02 cos(</)i -</)2) + VACos'^ff^^V'^^ 8A^] + SAp- 

where 8A^ is the reduction of the fundamental, 
given by Equation (22). 

* The theory contained in this section is similar to that 
given by King'^ but differs in detail because King's 
expressions are based on a cross-modulation formula of 
Huxley and RatcJiffe' which neglects the reduction of 
the carrier amplitude of the wanted wave caused by the 
di sturbing wave. 



18.8 dB 
(23) 



wave would do separately. If the arguments con- 
tained in the previous sections are followed on 
this basis we arrive at the result 



8A{t) =8Aum +8A^yjf{tf 



(25) 
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where 8A is the total increase in attenuation 

dA\jf(t) is the increase caused by the 
wanted wave on its own 
^^uf(0' is the increase caused by the un- 
wanted wave on its own* 

The received signal V{t) is again given by 
Equation (21), but 8A{t) now has the value given 
by Equation (25). Expansion of Equation (21) then 
shows that the fundamental component of the wan- 
ted modulation is reduced by 8Af + Mp, given by 
Equations (18) and (22). The harmonic distortion 
levels of the wanted transmission are given by 
Equations (23) and (24) but with 8Af replaced by 
8Af +MV. 

The modulation of the unwanted transmission 
is superimposed on the wanted transmission, to- 
gether with its second harmonic. The cross-modula- 
tion levels derived from Equation (21) are: 

Fundamental 



l'i = 20 logio(M' COS01 8Au) +8A^ + sa^^ - 12-8 



dB 



(26) 



Second harmonic 



Ls -20 logio 



((M')'cos02 5A'u) 



+ 5Ap+Mp-24.8dB 
(27) 

These levels express- the absolute modulation of 
the wanted transmission caused by the unwanted 
transmission; for example, a level of -20 dB would 
correspond to 10% modulation of the wanted trans- 
mission by the unwanted transmission. The levels 
also depend to a small extent on the depth of 
modulation of the wanted transmission, but this is 
a second-order effect; Equations (26) and (27) 
give the levels occurring when the wanted trans- 
mission is unmodulated 

Intermodulation frequencies f^ ± /^ and f^ ± 
2fJn, where f^ and f^ are the modulation frequencies 
of the wanted and unwanted transmissions respec- 
tively, are also generated. The intermodulation 
levels given by Equation (21), relative to the level 
of the wanted modulation, are: 

Frequencies f^ ± f'^ 

Lii=20logio(M'cos^i My) +8A^ + 8A^^ -18.8dB 

(28) 

Frequencies f^ ± ^fm 
Li2 = 2Ologio((M')^cos02 5A'u ) + 5A p + 5^^-30 .8 dB 



({MYcosc/y^ 5A'uj 



(29) 



Radio-frequency intermodulation products, such 
as the sum and difference of the two carrier fre- 
quencies, are not generated because the ionospheric 
time constant 1/Gvo is very large compared with 
the period of one cycle at r.f. No trace of these 
frequencies was detected in the recent experiments 
with high-power European transmitters^^ although 
it was expected that they would be observed. 



Throughout this section primed symbols refer to the 
unwanted wave. 



7. CONCLUSIONS 

When a powerful disturbing wave traverses the 
ionosphere the energy balance of the free electrons 
is upset and their collision frequency increased. 
As a result the attenuation suffered by this and 
other waves traversing the same region of the ionos- 
phere may be increased. The attenuation increase 
is proportional to the power of the disturbing trans- 
mitter and inversely proportional to its frequency 
unless it excites the extraordinary wave strongly, 
in which case the attenuation is accentuated near 
the gyromagnetic frequency. 

If the disturbing wave is amplitude modulated 
the attenuation will vary over the modulation cycle 
and the modulation will be impressed on all waves 
traversing regions of the ionosphere where signifi- 
cant attenuation changes occur. Medium-frequency 
waves are modulated in this manner in the E-layer 
but not in the D-layer. The cross-modulation is 
greatest at low audio frequencies because the 
changes in collision frequency do not occur instant- 
ly but are subject to a time constant of about 1-25 
milliseconds. 

If the disturbing wave itself is being received 
the attenuation increase will result in a non-linear 
relationship between the transmitted and received 
powers and the modulation waveform will be dis- 
torted. Ionospheric non-linearity may also give 
rise to audio-frequency intermodulation products 
if two powerful transmitters illuminate a common 
region of the ionosphere but radio-frequency inter- 
modulation products will not be generated because 
of the large ionospheric time constant. 

In view of the importance of these effects 
where high-power transmitters are concerned, expres- 
sions have been derived and are given in this 
report for: 

(a) the reduction and harmonic distortion of the 
received signal, caused by non-linearity 

(b) the cross-modulation suffered by all waves 
traversing the region of the ionosphere illumina- 
ted by a high-power transmitter 

(c) the cross-modulation and intermodulation which 
arises when two high-power transmitters are 
situated in the same locality. 
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